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Abstract 

We perform a rigorous mathematical analysis of the bending modes of a linear triatomic 
molecule that exhibits the Renner- Teller effect. Assuming the potentials are smooth, we prove 
that the wave functions and energy levels have asymptotic expansions in powers of e, where 
is the ratio of an electron mass to the mass of a nucleus. To prove the validity of the expan- 
sion, we must prove various properties of the leading order equations and their solutions. The 
leading order eigenvalue problem is analyzed in terms of a parameter b, which is equivalent to 
the parameter originally used by Renner. For < 6 < 1, we prove self-adjointness of the leading 
order Hamiltonian, that it has purely discrete spectrum, and that its eigenfunctions and their 
derivatives decay exponentially. Perturbation theory and finite difference calculations suggest 
that the ground bending vibrational state is involved in a level crossing near b = 0.925. We also 
discuss the degeneracy of the eigenvalues. Because of the crossing, the ground state is degenerate 
for < 6 < 0.925 and non-degenerate for 0.925 < 6 < 1. 



*This research was supported in part by National Science Foundation Grant DMS-0600944 while at Virginia Poly- 
technic Institute and State University, and also by the hrstitute for Mathematics and its Applications at the University 
of Minnesota, with funds provided by the National Science Foundation. 

1 



1 Introduction and Background 



In their original paper, [1], Born and Oppenheimer let be the ratio of the electron mass to 
the nuclear mass and expanded the wave functions and eigenvalues of the time independent 
Schrodinger equation in powers of e. We shall refer to such an expansion as a Born- Oppenheimer 
expansion. Since e is small, the first few orders of the expansions are thought to provide reasonably 
accurate results for the bound states of the molecular system. Often only the lowest (or leading) 
order terms of the expansions are even considered. 

The focus of this paper is the Renner- Teller effect (also called the Renner effect), which is later 
described in more detail. In short, a symmetry induced degeneracy exists in the electron states 
at a particular nuclear configuration, but when the nuclei move away from this configuration the 
degeneracy splits. As a result one must use more than one electronic state when attempting 
to solve for the total wave function and energy using the Born-Oppenheimer approximation. 
This effect was first predicted in 1933 by Herzberg and Teller [10] and was analyzed one year 
later by Renner [20] in a simplified model. We consider the current paper as an extension of 
the mathematically rigorous works related to the Born-Oppenheimer approximation, such as 
[3, 4, 5, 6, 7, 8, 13], to the model originally considered by Renner [20]. The main results are 
contained in theorem 2.1. We show rigorously that a Born-Oppenheimer expansion exists to all 
orders of e, with minimal mathematical assumptions. We prove that under our hypotheses, the 
molecular energy and wave function can be approximated by an asymptotic series in e that is 
truncated at arbitrary order. The leading order equations we obtain are unitarily equivalent to 
those found by Renner in [20]. This is the first rigorous derivation of the leading order equations 
of which we are aware. We feel it is especially important to make contact with a rigorous Born- 
Oppenheimer expansion here, since the Renner- Teller effect is not a straightforward application 
of the Born-Oppenheimer approximation. In their extensive review of the subject [16], Peric and 
Peyerimhoff give several interpretations of the origin of the Renner- Teller effect, and in particular 
they state "from the quantum chemical standpoint, the R-T effect is a consequence of violation 
of validity of the Born-Oppenheimer approximation." We will see that in the Renner- Teller case 
there is a valid Born-Oppenheimer expansion, but it differs significantly from the usual Born- 
Oppenheimer approximation since the degeneracy cannot be ignored. It must be analyzed in 
terms of degenerate perturbation theory. 

In recent years there have been several mathematically rigorous results justifying the validity 
of Born-Oppenheimer expansions under various hypotheses. The first rigorous proof related to the 
Born-Oppenheimer approximation in a physically realistic model was given by Combes, Duclos, 
and Seller [3, 4]. They proved the validity of the fourth order approximation for the eigenvalue 
and the leading order approximation for the eigenfunction. A few years later, Hagedorn proved [5] 
the existence of a Born-Oppenheimer expansion to all orders using the method of multiple scales, 
assuming that the potentials are smooth functions. In particular, he proved that for arbitrary K, 



2 



there exist quasimode energies of the form Ek{^) = X^^o E^^'^ and quasimodes of the form 
^k{() = Ylik=Q e'^ ^e'^^ that asymptotically approximate an exact eigenvalue and eigenfunction 
below the essential spectrum of a Hamiltonian H{e), in the sense that 

\\H{e)-i>K{e) - EK{e)^K{e)\\ < Ck e''+\ (1.1) 

The first five orders of E{e) were determined explicitly, and it is discussed how one could 
proceed to any arbitrary order K. These results were then extended to the case of Coulomb 
potentials for diatomic molecules in [6] and to general polyatomic molecules by Klein et al. [13]. 
Here, we will assume that the potentials are smooth, but we believe our results can be extended 
in a similar manner to the case of Coulomb potentials. 

2 Description of the Model and Statement of the Main 
Theorem 

Consider a triatomic molecule and fix the reference frame so that when the molecule is in the 
linear configuration, the middle nucleus is at the origin and the z-axis passes through all three 
nuclei. Let (0,0, and (0,0, -R2) be the coordinates of the upper and lower nuclei (so i?i > 
and i?2 < 0). We consider the bending modes by clamping the upper and lower nuclei to their 
fixed positions on the z-axis and allowing the middle nucleus to move in the perpendicular plane. 
Let (x,y,0) be the cartesian coordinates of this middle nucleus, and let (/O, (^) be the usual polar 
coordinates associated with {x^y) (see Figure 1). 
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Figure 1: The reference frame for the middle nucleus. 

If (xi,X2,--- ,X]\f) are the N three-dimensional electron coordinates, the electronic Hamilto- 
nian is 

1 ^ 

y) = - 2^ ^x, + V{x, y;xi,X2,--- , xn), 
i=i 
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where we have taken the electron mass to be 1, and the potential V includes the repulsion forces 
between the nuclei, the attraction forces between the nuclei and electrons, and the repulsion forces 
between the electrons. We think of h{x,y) as having parametric dependence on (i.e. it is a 

mapping from to the linear operators on the electronic Hilbert space), and we assume it is a 
real symmetric operator. We assume that y is a smooth function in all variables. Let be the 
ratio of the mass of an electron to the mass of the middle nucleus. Then, the full hamiltonian of 
this model is given by 

H{e) = - - A^,,y + h{x,y). (2.2) 

Let Tinuc = L^{M^, dxdy) and 'Hg? = L'^{1R^^), so that H{e) acts on the Hilbert space 
"Hnuc ® "Hei- We denote the inner product and norm on TCei by ( •, and 1 1 - j |g; and similarly on 
Ti-nuc by ( •, • )nuc and |H|„„c. 

Let Lf and L^"'^ = — ^ ^ be the operators associated with the projections of the electronic and 
nuclear angular momenta on the z-axis, respectively. The operator of total angular momentum 
about the z-axis is denoted by L™'^ = (/ Lf) + (L™'^ (g> I). We note that H{e) commutes 
with L™-^. We consider the electronic states when {x,y) = (0,0). In this case the electronic 
hamiltonian h{0, 0) commutes with Lf since the nuclei are in a linear arrangement. So, for \lf\ ^ 
there are two-fold degenerate electronic vectors ipi, Tp2 ^ Ti-ei satisfying /i(0, 0) -01 = Eq "01 arid 
/i(0, 0) V'2 = -E'o^ ' ^2) where Lf tpi = if ipi and Lf Tp2 = —if 1^2- Then, if the molecule is bent 
so that (x, y) ^ (0, 0), this degeneracy splits since the nuclei are no longer in a linear arrangement, 
and h{x, y) no longer commutes with Lf (see [12] for a discussion directly relating the breaking of 
symmetry with the breaking of the degeneracy). This is the Renner- Teller effect. As previously 
mentioned, the application of the Born-Oppenheimer approximation is not straightforward in 
this case. There have been numerous papers related to the Renner- Teller effect, few of which are 
relevant to our analysis here. We highlight one such paper by Brown and J0rgensen [2] for its 
completeness, and because it does discuss effects beyond the leading order. We encourage the 
reader interested to learn the historical development and recent findings of the theory to consult 
the review by Peric and Peyerimhoff [16]. 

Note that since changes in (j) correspond to an overall molecular rotation, the eigenvalues of 
/i(x, y) are independent of (j). Corresponding to the situation above where the electronic states 
at p = are linear combinations of eigenstates of Lf with eigenvalues if , —if 7^ 0, consider a 
pair of electronic eigenvalues Ei{p) and E2{p) of h{x,y) that are degenerate at p = 0, but the 
degeneracy breaks when p ^ 0. We refer to two such electronic states as an R-T pair with value 
\lf\. The eigenvalues of h{x,y) provide the usual potential energy surfaces for the nuclei, and 
there are several qualitatively different possibilities where the Renner Teller effect is important. 
See Figure 2. We refer to [12, 14] for further examples and discussion of Renner- Teller surfaces. 
We focus strictly on an R-T pair of states corresponding to sketch (a) of Figure 2, where both 
surfaces have local minima at p = 0. In this case the optimal nuclear configuration, corresponding 
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to both electronic states of the R-T pair, is hnear. This was the situation considered by Renner 
[20] in 1934. 




P 



Figure 2: Potential energy surfaces of three qualitatively different cases corresponding to an R-T pair 
of electronic states. 

Throughout this paper, we assume the following hypotheses: There is an R-T pair of states, 
with eigenvalues Ei and E2 both having minima at p = 0. We assume that for some neighborhood 
of p = 0, El and E2 are isolated from the rest of the spectrum and are in (x, y). This implies 
that Ei{p) and E2{p) have asymptotic expansions in powers of p^. We assume that splitting 
occurs at 2nd order, that is, that Ei and E2 are asymptotic to ^-^-^ p'^ and p^ for small p, 
respectively, for some <b < a (we have taken Ei{0) = -©2(0) = for convenience). Renner [20] 
argued that an R-T pair with value \lf\ = 1 will exhibit splitting at 2nd order, an R-T pair with 
value \lf\ =2 will exhibit splitting at 4th order, and in general an R-T pair with value \lf\ = n 
will exhibit splitting at order 2n. We instead assume 2nd order splitting occurs and later prove 
that the R-T pair has value \lf\ = 1, agreeing with Renner's argument. We are now ready to 
state our main theorem. 

Theorem 2.1. Assume the hypotheses described above, in particular that the potentials are smooth 
and there is an R-T pair Ei{p), E2{p) that are asymptotic to ^^p^ and ^^p^ respectively, where 
< b < a. Then for arbitrary K, there exist quasimode energies 

E,,K = ^^'^ 

k=0 

and quasimodes 

K 

that satisfy 

\\{H{e) - E,,k) < CKe^^' 1 1 1 • 

The quasimodes are associated with the local wells of Ei and E2 in a neighborhood of p = 0. 
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Remarks: 



1. Quasimode estimates correspond to discrete eigenvalues of H(€) when -E^./^ lies below the 
essential spectrum as characterized by the HVZ theorem [19]. 

2. Since [L™-^,i?(e)] = 0, quasimodes can be constructed to be eigenfunctions of L™-^, 
with eigenvalues Z™"^ G Z. The eigenstates of H{e) corresponding to Z™"^ = are non- 
degenerate, while the eigenstates corresponding to |Z™"^| 7^ are two- fold degenerate. In 
particular, there is an Z™"^ state and a — Z™"^ state, with eigenfunctions that are complex 
conjugates of one another, together forming a degenerate pair of states associated with H{e). 

3. The first two orders E'^'^'^ and E^^'^ are zero and the second order E^"^^ is determined by the 
leading order eigenvalue equation = E^'^^'^, on the Hilbert space 



4. The presence of the two levels Ei{p) and E2{p) gives rise to twice the number of vibrational 
levels as usual in the following sense: If 6 = 0, the upper and lower component equations 
of the leading order equation which determines E^'^\ are both two-dimensional harmonic 
oscillator equations. So, there will be two eigenfunctions, one associated with the each of 
the upper and lower components, for each of the usual eigenstates of the usual harmonic 
oscillator. Then for small 6, this will give rise to two vibrational states via a perturbative 
approach, for each of the usual harmonic oscillator states. This is shown in detail in section 



5. For small h and e, the ground state of H[e) (meaning the lowest vibrational level correspond- 
ing to the R-T pair we are considering) is degenerate, corresponding to a pair of states which 
are eigenfunctions of L^*-^^ with eigenvalues t^^'^ = ±1. In section 5, we give plots which 
suggest that for approximately 0.925a < b < a the ground state is non-degenerate, corre- 
sponding to a state with Z™-^ = 0. 

The paper is organized as follows: In section 3, we derive perturbation formulas to construct 
the quasimodes that will enter in our main theorem. In section 4, we prove various properties 
of the leading order Hamiltonian that are needed to prove the main theorem. In section 5, we 
analyze the leading order eigenvalue problem. Only some of the eigenvalues and eigenfunctions 
of the leading order equation are solved for exactly. In section 6, the degeneracy structure of the 
full Hamiltonian H{e) is discussed. In section 7, we use the results of the previous chapters to 
prove the main theorem. 



L2(iR2^ dXdY;<C?), where 




The higher order E^^'^ are determined through the perturbation formulas presented in chap- 
ter 3. All odd order E^^'^ are zero (see Appendix). 



5. 
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3 The Construction of the Quasimodes 



Before we begin the formal expansion, we first fook at some properties of the electronic eigenvec- 
tors and eigenvalues, construct electronic basis vectors that are smooth in terms of the nuclear 
coordinates, and derive the leading orders of the matrix elements of the electronic Hamiltonian 
in this basis. 

3.1 The Two-Dimensional Electronic Basis Vectors 

For the N electrons, as well as the nuclei, we use the same fixed reference frame previously 
described. Let {rj, 9j, zj) be the cylindrical coordinates of the jth electron in this frame. Suppose 
that for p > 0, ip{x, y; 6i, 92, - ■ ■ , 9n) '■ JR^ T~iei is an electronic eigenvector of /i(x, y). We have 
suppressed the dependence on rj and Zj because it is irrelevant to the discussion here. The 
electronic eigenfunctions are invariant with respect to a rotation of the entire molecule. So, the 
eigenfunctions have the property 



for p > 0. It follows that if ij){x,y) is continuous at p = 0, then -0(0,0) has no 9j dependence. 
Since an eigenvector corresponding to an R-T pair with positive value will have some 9j 
dependence at p = 0, we do not have well-defined continuous electronic eigenfunctions of y) 
in a neighborhood of p = 0, that correspond to an R-T pair with value |?^'| > 0. We need basis 
vectors for the two-dimensional eigenspace of Ei{p) and E2{p) that are smooth in x and y. The 
matrix elements of h{x^ y) in our electronic basis determine the form of the leading order equations 
to follow. We note that in deriving these matrix elements, we do not use the matrix elements 
of Only second order splitting in Ei and E2 is needed, as well as the fact that our smooth 
basis vectors are not eigenvectors of h{x,y). This gives rise to off-diagonal terms in the basis 
representation of h{x,y). In this sense, the unusual form of the leading order equations can be 
thought of as a result of the discontinuity of the electronic eigenvectors in the nuclear coordinates, 
i.e. there is no smooth electronic basis that diagonalizes the electronic hamiltonian. We note that 
matrix elements we derive here, are related by an {x, y)-independent unitary transformation to 
those given by Yarkony [22]. See also Worth and Cederbaum [21] for a general discussion of the 
topology and classification of different types of intersections of potential surfaces. 

We now describe our approach. Choose any two normalized orthogonal electronic vectors ipi 
and 'ip2 that span the eigenvalue eigenspace of /i(0, 0). Let P{x, y) denote the two dimensional 
projection onto the electronic eigenspace associated to the two eigenvalues of h{x, y). For small 
X and y, define 



i){x, y; 6*1, 6*2, • • • , 



9n) = ^{p, 0; 9i-^,92-<P,--- , 



On-^) 



y) 



1 



P{x, y)t(;i. 



(3.1) 
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Let Pi{x, y) denote the orthogonal projection onto this vector, i.e., 

Pi{x, y) = \^i{x, y)){^i{x, y)\. 



Next, define 



and 



x{x, y) = (1 - Pi{x, y)) P{x, y) ip2, 



'^2{x,y) = ==xix,y). (3.2) 

V \ x{x, y), x{x, y) ) 

Then { 4'i(a;, y), ^2{x, y) } is an orthonormal basis for the range of P{x^ y). From the formula 
[19] ^ 

P{x,y) = j^[\-h{x,y))~U\ 

where C is a closed path in the complex plane encircling Ei{p) and E2{p) but no other spectrum 
of h{x,y), we see that these vectors are smooth in x and y, since we have assumed that the 
potentials are smooth and hence the resolvent of h{x,y) is as well (recall we are only working in 
a neighborhood of the origin {x,y) = (0,0)). Note that we can arrange for these vectors to be 
real, which we assume has been done. 



3.2 The Matrix Elements of the Electronic Hamiltonian 

The span of { ^'i(x, y), ^2(2;, y) } is an invariant subspace for h{x, y). Using coordinates in this 
basis, the restriction of h{x, y) to this subspace is unitarily equivalent to the real symmetric 
matrix 

hu{x, y) hi2{x, y) 

h2i{x, y) h22{x, y) 
where 

hjk{x,y) = {^j{x,y),h{x,y)^k{x,y)). 

Again, since we have smooth potentials, hij{x,y) can be expanded in powers of x and y. Since 
we assume the degeneracy splits at second order, the eigenvalues of this matrix are 
Ei{x, y) = {x'^+y'^) + 0{p'^) and E2{x, y) = ^ {x'^+y'^) + O(p^). Using these expressions 
for the eigenvalues we show that up to an (x, y)-independent unitary transformation, this matrix 
is 



a + b 2 o — 6n 

X + y ±bxy 

2 2 

a — b n a + 62 

±bxy X + y 

^ 2 2 

To show this, we consider a traceless, real symmetric matrix 

/iii(x, y) hi2{x, y) 
h2i{x, y) -hii{x, y) 



(3.3) 



(3.4) 
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with eigenvalues E±{x, y) = ^ + 0{p^). The form in (3.3) wiU fohow from the analysis below. 
Using (3.4), we have the characteristic equation 



El + 0{f) = hl^ + hl^. (3.5) 

By expanding in powers of x and y and equating orders in the above equation, it can be easily 
shown that the constant and linear terms of hn and hi2 must vanish. We then write. 



hii{x, y) hi2{x, y) 
h2i{x, y) -hii{x, y) 



Ax^ + By^ + Cxy + 0{p-^), (3.6) 




where A, B, and C are traceless 2 by 2 matrices with constant entries. We can apply a constant 
unitary transformation to (3.6) that diagonalizes A, which we assume has been done. An obvious 
consequence of (3.5) and E±{x, y) = ± + 0{p^) is that if A is diagonal, it must be 



A 



We let 

(fell bi2 \ I Cii Ci2 

, and C = 
bl2 -hi J \ Cl2 -cii 

and use (3.5) with E±{x, y) = ±fP + 0{p^) to solve for bij and Cij by equating the powers of x 
and y. This gives us four equations, the first equation comes from the coefficients, the second 
comes from the x^y'^ coefficients, etc. 

/ : 1 = fen + fe?2 (3.7) 

x^y^ : 2 = 2bn + cli + c?2 (3.8) 

xy^: = 2(fencii + 612C12) (3.9) 

x^y: = 2cii (3.10) 

These equations have 3 solutions. Two of the solutions are 

(fell,fel2,Cii,Ci2) = (-1,0,0, ±2), 

which give 

hii{x,y) hi2{x,y) \ f x^ - y^ ±2xy 
h2i{x,y) -hii{x,y) J \ ±2xy -(x^-y^) 
These solutions give rise to (3.3). The only other solution of equations (3.7)-(3.10) is 

(fell,fel2,Cll,Cl2) = (1,0,0,0), 
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which gives rise to 



hn{x, y) huix, y) 
h2i{x, y) h22{x, y) 



2 ^ 



\ 







We do not consider this case. Aside from being uninteresting, it imphes that the basis vectors are 
the eigenfunctions of y) (at least to leading order). We assume that the off diagonal terms 
in (3.3) are hxy, since the —bxy case is related by the trivial change of coordinates y ^ —y. 



3.3 The Formal Expansion 

To construct the quasimodes in theorem 2.1, we introduce the scaled variables {X, Y) = {x/e,y/e). 
The intuition of the Born-Oppenheimer approximation suggests that the adiabatic effects will 
occur on the {x,y) = {eX, eV) scale, whereas the semi-classical motion of the nuclei is determined 
on the {X, Y) scale. In terms of the {X, Y) variables, the Hamiltonian in (2.2) is 

H{e) = -jAx,Y + Hex, eY). 

We define 7i to be the Hilbert space L^{1R^, dXdY- €?) and we denote the inner product on this 
space by (•, 

We seek solutions to H{e) *&(e, X, Y) = E{e) ^(e, X, Y). The wave function ^'(e, X, Y) can 
be written in terms of the orthonormal basis functions { ^1(2;, y), ^'2(2;, y) } from (3.1) and (3.2) 
as 

^{e,X,Y) = f{e,X,Y) ^i{eX, eY) + g{e, X,Y) ^2{eX, eY) + i^^{e, X,Y), (3.11) 

where ( )ei = 0. 

Substituting (3.11) in H[e) ^'(e, X, Y) = E{e) ^'(e, X, Y) gives three equations; one along 
^1, one along ^'2, and one in span{'ifi, ^2}^- We denote the projection on span{'^i, 'I'2}"'" by 

P±. 
Along 

- — Ax,Y f + hu f + hi2 g - — { ^l, Ax,Y ^±)el - f { ^l,A^.,y ^1 )el 

/ T A TV 3 f dg , d^2 , dg ^ 8^2 , \ 

= E{e) f. (3.12) 

Above we have used that ( ^'j, )ei = 0, which we know from normalization and the fact 
that the electronic basis vectors were chosen real. Along \I'2 we get a similar equation with 
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f ^ g, ^1 ^ ^2, hii ^ /l22, hi2 ^ /l21- 

In spaniel, ^2}^: 

2 44 

- |- [ Ax,Y + {h P±) V± - J f P±[ ^x,y *i ] - Y [ ^-.s/ *2 



9/ 



dx 



df 



5y 



do 

+ — Pi 



d^2 

dx 



dX 
= E{e) V'^. 

We adopt the following notation for simplicity: 

Tij{x,y) = {^i,A^^y^j)eU 

y) 



dy 



(3.13) 



9x 



dy 



We have identities involving these quantities since {^i, ^2} are orthonormal and real valued. 
For instance we know the diagonal elements of A and B are zero and A12 = —A21, B12 = —B21. 
Now we expand all functions and operators with e dependence. For example, /(e, X, Y) = 

00 

^ f^^\X, Y) . For functions and operators with exclusively (x, y) dependence, we know the 



k=0 



form of the expansions. For example, ^1(2;, y) = ^i(eX, eF) = ^ e'' Y), where 



k=0 



^['\X, Y) = ^2 J(j^, d!i3yU (O'O) Equations (3.12) and (3.13) become: 



oc / ^ \ CO k 

fc=2 ^ fc=0 i=o 



"■12 y 



+ E ^' E 

fc=2 3=2 



{^^r^\ Ax,Yi^t'^ )el 



+ E^^E 

k=A j=4 
00 fc 

+ E^'E 

A:=3 i=3 



(i-3) 
12 



d 

dX 



B 



0-3) 
12 



d_ 
dY 



,(fc-j) 



^ ^i^O) /(fc-i) (3.14) 
fc=o i=o 
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and 



^ P 



oo k 
k=2 3=2 



oo k j 

A:=4 i=4 Z=4 



{i-2) 



oo k j 

fc=3 i=3 1=3 



Ax,Y 



(fc-i) 



oo k 
fc=0 j=0 



2 ■ • ^' 



-i) 



oo k j / 

+ E^'EE h^-- 

fc=3 j=3 Z=3 \ 



+ P 



dx 



(1-4) 

d 

dX 



P 



(i-0 



9y 



d_ 
dY 



P 



U-i) 



dx 



d 

dX 



P 



U-i) 



dy 



d_ 
dY 



(3.15) 



We now collect terms at each order of e. Recall there is an equation along ^'2 analogous to 
(3.14). At each order, we will combine these two similar equations into one matrix equation. 
Order The e'' terms require 



'ni "'12 

4? 4? 



,(0) 



fiO) 



(3.16) 



(3.17) 



The hij{x, y) vanish until second order, so this forces E^^^ = in (3.16), and consequently = 
after applying the reduced resolvent [( /i(x, y) Pi(x, y))*-°^ in (3.17). 
Order 1 As above, the terms reduce to 



/(O) 
n(0) 



0, 



12 



So we get = and = 0. 

Order 2 Using the known second order terms for the hij{x,y), the terms require 

/(o) \ / f{0) 



Ma 



{hPi 



)(0) ^(2) 



0, 



where 



V 



a + 6 



2 

6xy 



1 A a- 



6xy 
6 



9 a + b n 

X^ + ^^r^ 



Recall we have assumed the +hxy case for the off diagonal entries. By again applying the 

(2) 

reduced resolvent in the last equation we have = 0. In chapter 4 we show that EI2 is self- 
adjoint (on the correct domain) and has purely discrete spectrum with infinitely many eigenvalues 
for a > 6 > 0. We are only able to solve for some of them exactly. In chapter 5 we show that 
there is at most a two- fold degeneracy in the eigenstates of IHI2, but that no splitting occurs 
in the quasimode eigenvalues, i.e., the degeneracy remains to all orders of e. We can therefore 
proceed as if the eigenstates of BI2 were non-degenerate, since we can take any linear combination 
of degenerate states for f^'^^ and g^'^^ and we know it will lead to a valid quasimode and energy 
E{e). Fix E^'^\ f^'^^ and g^^^ corresponding to one of the states of Ha. 
Order 3 The terms require 

/(O) ^ ^ I f{i) I' f{0) 

(3.18) 



,(0) 



+ 



£;(2) 



^(3) 



.(0) 



p 



(0) 



dx 



(0) 



dX + ^ 



dy 



(0) 




where 
H3 



(0) 



+ P 



'ni "'12 

"'21 "22 



dx 



+ 



dX + ^ 



dy 



A. 



(0) 

21 



d 

dX 



dY 



d 



'2\ 



d_ 
'dY 



A 



(0) 



d 



B 



(0) 



d 



(3.19) 



12 dX 12 QY 




Since EI2 is self-adjoint, we can take inner products of both sides in (3.18) with 



/(O) 
n(0) 



obtain 



Ei^) 




to 
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In the appendix we argue that ah of the odd order E^'^^ are zero. Let Qj_ be the projection 
in H onto the subspace perpendicular to the eigenspace of the eigenvalue E^"^^ of EI2. Adopting 
"intermediate normalization" we may choose the non-zero order wave functions perpendicular to 
the eigenspace of E^"^^ (note that this will produce a non-normalized quasimode), so that 



f 



(k) 



for k > 1. Then from (3.18) we get 



9 



(1) 



From (3.19) we have 



(3) 



ihP±) 



(0) 



p 



(0) 



dx 



(0) 



dX 



+ P 



(0) 



dy 



(0) 



5_ 
W 



(3.20) 



f{0) 



+ P 



(0) 



dx 



d 

dX 



Order 4 The terms require 

;(2) 



9 



(2) 



+ 



i^(3) 



+ P 



(0) 



dy 



d_ 
dY 



+ 



f(0) 
,(0) 



(3.21) 



0, (3.22) 



{hP^f> = - {hP^f^ V? 



(A,,, M/i)(°) 1 + Pf [(A,,,vt2 



,(0) 



(0) 



4 i 

i=3 «=3 

+ E E ( 

i=3 Z=3 



where 



(i-0 



5^1 



(i-3) 



5x 
dx 



7^(0) ^(0) 

7^(0) rpiO) 



(1) 5 



9 



9*1 



9 



+ P 



(i-0 



+ 



"■11 "-12 



41' 4? 



dy 

d^\^'-'^ 
dy 

+ 



d_ 
'dY 



d_ 
dY 



,(1) 9 



"■12 



dX 



B 



(1) 5 
12 5F 



(1) 9 



21 



dY 



(3.23) 



;(4-.) 
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Using what we know through order 3, we can solve (3.22) and (3.23). From (3.22) we obtain: 



;(0) 



5« 



+ 



f(0) 
(0) 



n 



9 



H 



and 



From (3.23) we get 



(4) 





IH2- 




get 








-1 






r 



Q. 



i^(3) 



+ 



;(0) 
n(0) 



(A.,, /W + Pf (A.,, ^2) 



,(0) 



4 i 

+ E E I ^- 

i=3 Z=3 



0-0 



4 j 



+ E E I 

j=3 /=3 



^ , pU-i) 



5y 



(/-3) 



dY 



d_ 
dY 



i(4-j) 



Order > 5 We now show that we can proceed in this manner to any order of e desired. In 
chapter 4 we wih show that all of the quantities involved exist in the relevant Hilbert space. If 
k > 5, the terms require 



(fc-2) 



+ (m-s-E^^^ 



9 



+ E 

j=4 



fiO) 



k-3 

+ E 

i=2 



1^ / {^t"\ Ax,y Vr^'))e. 



0, 



(3.24) 
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fc-3 



Ax,Y ■0 



i=2 



fc-3 



j=4 /=4 



1 ^ "0-0 r^A ,T,_^a"4) 



k j 

+ E E I 

j=3 /=3 
k j 

+ EE 1^"- 

j=3 /=3 \ 
J=2 



(i-0 



0-0 



9x 



(1-3) 



5y 



(Z-3) 



9_ 

9y 



<-(fc-i) 



dx 



d 

dX 



P 



0-0 



5y 



9y 



where 



tO-4) ^0-4) 

7^0-4) ^0-4) 
-^21 ^22 









( /.Si' 















A 



. _ ,0-3) d_ _ „0-3) 5_ 

V ^21 5^ ^21 5^ 



+ 

0-3) d_ „(j-3) d_ 

12 dX 





for j > 4. 

Following what we have seen through order 4, assume from previous orders that 



/O) 
,0) 



for j = 0, 1, . . . , - 3, E^^^ and V?^ for j = 0, 1, . . . , A; - 1 



(3.25) 



are already determined. Then, we can solve (3.24) and (3.25) for Z^'^-^), c/C^-^), ^^^J"), and E^''\ 
From (3.24) we obtain: 



fc-i 
i=3 



;(0) 

T fc-3 

- y 

2 ^ 

J=2 



£;0) 



n{0) 



+ 



/(O) 
,(0) 



n 



(^(^-^\ Ax,Y^^t'^ )el 



(3.26) 
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and 



fiO) 



k-1 



f{k-j) 



1 ^3 / ( Ax,Y hi 



(0) 



From (3.25) we get 



(k) 



ihP±) 



(0) 



g 



k-3 



, , ^,0-2) A.... jS^-j) 

- y 

^ U \ { ^^^^ Ax,y 



(3.27) 



k-3 



J=2 



(i) 



+ E E ^ ( ^i)'"'^ 1 f''-'^ + ^^"'^ 



i=4 «=4 

+ E E ( ^- 

j=3 1=3 

k j 

+ E E I 

j=3 1=3 



U-i) 



dx 
9x 



5y 



^ ] fik-j) 



<9y 



9_ 

9y 



+ E E^'^ i^V^ 
i=2 

So we can proceed in this manner to obtain ^'(e) and E{e) up to any order in e. 



-,{k-j) 



(3.28) 



4 Properties of the Leading Order Hamiltonian 



We adopt the foUowing notation throughout: 
/ 1 ' 



1. We let h 



1 



. If yl is an operator on the Hilbert space L'^{IB?, dXdY), then 



A /2) is the operator on Ti. given by 



A 
A 



2. If D{A) is the domain of the operator A on the Hilbert space L^{1R^, dXdY), then 
D{A ® h) = D{A) e D{A) C n. 

In what follows, we prove various needed properties for the expansion to all orders. Let 



12 = - - /\x,Y ® h + 



( a + b ^2 , « - ^ 
2 2 

bXY 



bXY 

a-b ^2 ^ a + b ^2 
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Note that if we let (X, Y) = (a^/^ X, a^/^ and 6 = - , then 

a 



+ b 



2 

V 



X^ - 



XY 



XY 



- [X^ -Y^ 



We now use the Kato-Relhch Theorem [18] to prove self-adjointness of IHI2. 



(4.1) 



Theorem 4.1. If a > b > 0, then M2 is self-adjoint on Dho © Dho j where Dho is the usual 
Harmonic oscillator domain in L"^ (IR^ , dX dY) , and essentially self-adjoint on Dho © Dho ; 
where Dho is o-ny core for the usual Harmonic oscillator. 



Proof: 

Define 

and 



Hho 



V{b) 



- (X^-Y^] 



XY 



XY 



1 



We prove that for < 6 < 1, V{b) is relatively bounded with respect to Hho, with relative bound 
b. The conclusion then follows from the Kato-Rellich theorem [18] and (4.1). 

For each fixed X and Y, the eigenvalues of V{b) are it | {X"^ + Y"^). It follows that 



V{b)v 



< b 



- {X^ + Y^)® h]v 



where v € is any two component vector, and we use the usual Euclidean norm. This inequality 
implies the L^(M^, dX dY; C^) = H norm estimate 

1 



y(6)V 



n 



< b 



{X^ + Y' 



H 



where ip{X, Y) & TC is a two-component vector-valued function. 
We now show that 



V(b) i) 



n 



< b 



i (X^ + Y^)0l2]iJ 



< b II HHoiPWn + ^ II "0 llw- 



(4.2) 



n 



for all ijj G Dho © Dho ■ We have already shown the first inequality. The hard part is the second 
estimate, which follows from 



{{X^ + Y^) (g) l2)ip\\ < II ((-Ax,y + ^2 + ^ j^) Vll + 2 II VI 
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This easily follows from 

+ ® Is) < \\[[-Ax,Y + X^ + Y^)^l2)'ipf + 4\\i;f. (4.3) 

Rather than proving this directly, let us first prove a simpler relative bound estimate for the 
operators on L^{1R, dx). We show that for (j) ^ D 



92 



+ X' 



dx"^ 



+ x' 



+ 2 



(4.4) 



To prove this, let p = — i^, and calculate the commutators 



[x, p\ 



and 



[X, p 



2ip. 



2 



< 



(4.5) 



We have 

{<P,xU) 

{cP, ((/ + x2)2-xV-pV -/),/,) 

||(/ + 2;2) (f)\f - {(p, (2;V+/a;2) 

In this last expression, we use the commutators above to write 

{(p, [x^p^ + p^x"^) (p) = {4>, [xp^x + p^] + xp^x + [p^, x]x) cp) 

= 2 {(p, xp^x (p) + 2i {(p, {xp — px) (p ) 

= 2{(p, xp^xcp) -2 {(p, (p). 

In this last expression, the first inner product is the expectation of a positive operator (since xp^x 
has the form A* A with A = px). Using this and (4.5), we see that 

\\x^(pf < Wiv' + x^) (Pf + 2 \\(P\\\ 

and (4.4) is proved. 

Now we simply mimic the proof of (4.4) to prove (4.3). We write 

II (X^ + y"^) (/.f 

= (0, ((-Ax,y + X2 + y2)2 _ ^2^^^ ^ ^^^^ ^^2 ^ y2^ ^ (^2 ^ y2^ Ax,y)</>) 

The operator y is positive. The operator — Ax Y'^ = —Y^ Ax is also positive since it equals 
A* A with A = pxY- Similarly, — Ay X"^ = — X"^ Ay is positive. By the commutator tricks we 
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used above, — Ax X'^ — X"^ Ax and — Ay — Y"^ Ay each are positive operators minus twice 
the identity. Thus for all (j) € Dho-, 



Ax,Y + X^ + Y^f - Aj^Y + ^x,Y {X'^ + Y'^) + {X^ + Y'^) Ax, 



Y 



< 



and hence, 

||(X2 + y2)^||2 < ||(_Ax,y + X2 + y2)^||2 ^ 4 || || ^ . 

It follows that (4.3) holds for all G Dho © Dho ■ This proves (4.2) and the theorem follows. □ 



Unless otherwise stated, it is assumed that by H2 we are referring to this operator with domain 
D(J\l2) = Dho © Dho- We now show that IHI2 has purely discrete spectrum. 



Theorem 4.2. If a > b > 0, EI2 has purely discrete spectrum, with countably many eigenvalues 
{1^3(^2)} jLi satisfying 



NVa- b < /iAr(Ar„i)+i(E[2) < /iAr(Ar-l)+2(El2) < ... < ^Ar(Ar+l) (IHI2) < Ny/a + b, 

for N = 1,2,3,... 
Proof: 

Let (p, (j)) be the usual polar coordinates associated with {X, Y). Define the unitary operators 
U, W : 7i ^ TC hy (defined as multiplication operators on 7i): 



U 



cos((/)) — sm{(j)) 
sin((/)) cos((/)) 



and W 



1 

V2 



Define 



U"^ Mo U 



and 



/ 1 



/ 1 . 1 a+b 2 

_ ]_ d_ 

d(f) 



1 a lb 6 9 



2 "''^ ' 2p2 ■ 2 
d(f) 



P 



1_ d_ 

,2 PiA 



\ 



1 a-6 



1 a±6 



2 



1 ^ 

2 ^'-^ ^ 2^92 + 2 



and note that Hq < Hq < HfJ". Now we define 
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In the context of the min/max principle [19], for all n G N, 



/ 1a a ± 6 n 







2 ^P.'/' + 



a±6 




2 



/i„(]H-) = ^„(]Ho) < fXnim) = fXnim) = ^n(IHIo) < /in(IH+) = ^in(IH+). 



The operators have purely discrete spectrum, with 2A^-fold degenerate eigenvalues of a ±b 
for N = 1,2,... So, EI2 must have purely discrete spectrum with eigenvalues /ii(E[2) < fj,2i^2) ^ 
. . . satisfying the required bound. □ 

To prove the quasimode can be expanded to any order in e, we must show the terms arising 
at arbitrary order in the equations of chapter 3 are in TC. This follows from the propositions 
and lemmas we now prove. A similar analysis was needed in [8] and the proofs presented here 
are analogous to those found in [8]. For our purposes it must be shown that the details can be 
extended to this situation on 7i. 

Lemma 4.3. Let T{a) be defined on a dense domain of a separable Hilbert space H, and suppose 
that T{a) is an analytic family in the sense of Kato for all a G C, and self-adjoint for all a £ M. 
If, for all a G M, T{a) has purely discrete spectrum with eigenvalues accumulating at 00, then 
T{a) has purely discrete spectrum for all a G C. 



First we note that if a self-adjoint operator has purely discrete spectrum with eigenvalues 
accumulating at 00, then it has compact resolvent by Theorem XIII. 64 of [19]. We also note that 
for any closed operator A, {A — is compact for some /.t G p{A) if and only if {A — is 
compact for all ^ G p{A). This follows from the first resolvent formula. 

We first show that if T{a) has compact resolvent for all a G Si, then T[a) has compact 
resolvent for all a G C. We then show that if a closed operator defined on a separable Hilbert 
space has compact resolvent, then it must have purely discrete spectrum. 

Since T{a) is an entire analytic family, the resolvent Ra{^) = (^(a) — A)~^ is analytic in 
both a and A inside the set R = {(a, A) : aGC, AG p{T{a)) }. From Theorem XII. 7 of [19], 
R is open in both a and A. Let 53(-ff) and C(-ff) denote the bounded operators and compact 
operators on the Hilbert space H respectively. It follows from the Hahn-Banach Theorem [17], 



that for any B G *B(if) \ <1{H), there exists Ib G «B(if)* such that Ib{B) / 0, and = on 



Proof: 



<t{H). 
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Note that since T{a) is an analytic family, we know the resolvent set is non-empty for all 
a G C. Define the set 

T = { a G C : -Ra(A) is compact for all A G p{T{a)) }. 

We show that T = C. 

Let Bs{z) denote an open disk in the complex plane of radius s > 0, centered at z G C. 
Let Aq G p(T(0)). Since the set R is open, we know that there exists a disk -65(0), such that 
Ao G p{T{a)) for ah a G ^^(O). Let / G *B(if)*, such that / is vanishing on ^{H). Then 
the function /(a) = l{RaiXo)) defines an analytic map from -65(0) into C. Since the resolvent 
is compact for a G JR, we know f{a) = for all —6 < a < 5, which implies /(a) = for all 
a G Bs{0)- Since I was chosen arbitrarily in { / G ^{H)* : / vanishes on <t{H) }, it follows that 
-Ra('^o) is compact for all a G Bs{0) (if Ra'{Xo) was not compact for some a' G -65(0), we could 
choose / so that (Aq)) 7^ 0, contradicting /(a') = 0). Hence .65(0) C T. 

We now assume that T 7^ C and show this leads to a contradiction. Let 
r = sup{(5 > : Bs{0) C T}. Note that < r < co since we have assumed T / C. Then, 
there exists with |ao| = r, such that every neighborhood of oq contains a point not in T. Let 
Ao G p(T{ao)) and choose 6' > small enough so that Aq G p{T{a)) for all a G Bs'iaQ). Choosing 
I as before, we know g{a) = l{Ra{Xo)) is analytic on Bsi{ao) and g{a) = on Br{0) n Bs' (ao). 
So, g{a) = on all of ^^/(ao)- Again since / was chosen arbitrarily, there exists an entire 
neighborhood of ao in T. This is a contradiction, so T = C. 

We now show that a closed operator with compact resolvent has purely discrete spectrum. Let 
A be a closed operator. Fix A G p{A) and let R{X) = {A — A)~^ be compact. Then the spectrum 
of R{X) is made up of at most countably many eigenvalues of finite multiplicity that can only 
accumulate at [17]. For E ^ X, we have 

A-E = A-X-{E-X) = {E- X){A - A) ( - R{X) 

\E — X 

From this we see that if ;gr^ G p{R{X)), then E G p{A). So, if ^ G a{A), then ^ G (t{R{X)) 
and thus -^^^ is an isolated eigenvalue of R{X) with finite multiplicity. Since, 

i?(A)^' = ^ <^ - A) ^' = - A) ^' <^ A"^ = E"^, 

E — X 

it follows that E is an isolated eigenvalue of A with finite multiplicity. Therefore, cr{A) is made up 
of at most countably many eigenvalues of finite multiplicity that can only accumulate at infinity. 
The conclusion of the Lemma follows. □ 

Before we prove Proposition 4.4, we consider a different decomposition of EI2. We define Hq 
and V to be 

/a + br,a — bn \ 



Ho = — - Ax,Y ^ h and V 



+ bXY 

2 2 

a — b o a + b o 
bXY X^ + —— Y^ 
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so that IHI2 = -ffo + ^- Note that for any X, Y, the eigenvalues of V are ^ {X^ + Y"^) and 



2 

(X2 + y2)_ So for f, g£ L^{]R^) 



and y is a positive operator. 



Proposition 4.4. Lei ^' = [ ) G W 6e a solution of = , with E > 0. Then, 

\9 J 

f, ge C°°{M^), V/, Vg G L'^{R'^), and for any 7 > 0, 

/, g G Z)(e^<-)), V/, V<7 G Z?(e^<->), A/, G Z)(e^<-)), 
w/iere (x) = \/l + + 
Proof: 

Let Vu = + ^y2^ ^12 = 1^21 = bXY, and F22 = ^X"^ + ^Y^. 

Then, /, g satisfy the following pair of equations: 

{-A + Vu) f + Vi2g = Ef (4.6) 
(-A + ^22) 5 + ^21/ = Eg (4.7) 

To show that /, 9 G C°°{1R^), we follow the proof of Theorem IX.26 of [18]. Let O be a bounded 
open set in M^. Since f,ge L'^{1R^) = Wq and the Vij G C°°, we have Vuf, V21/, ^125, ^^225 G 
Wo{n). It follows from (4.6) and (4.7) that A/, Ag G WoiQ). Then by the Lemma on pg. 52 
of [18], f,g£ W2{^)- Repeating the argument we get f,g£ Wm{^) V m G It follows from 
Sobolev's Lemma that f,g€z C°° on Q. Since O was arbitrary /, 5 G C°°(iR^). 

We now show V/, Vg G L^. We know ^' G 15(1112). Let D{-A) and Q(-A) be the domain of 
self-adjointness and quadratic form domain of —A respectively. Then 

D{M2) C D{-A) e D{-A) c Q(-A) © Q(-A) = | ^ = ^ ^ ; y/, Vg G L^{1R^) | . 

We now use the Combes-Thomas argument (see theorem XIII. 39 of [19]) to prove that 
f,gG D{e"'^^^). The argument can be repeated for D{e"'^'^^), and since 

e7M < e7 eTd^l + l^-l) < g7 e27max{|X|, < g7 j^g27|X| _^g27|y|^ ^ 

we then have f,gE D{e^^^^). 

For a G M, consider the unitary group W{a) = e*"'''' (8> I2 and the operator 
M2(q) = l^(a)M2T^(a)-i. We have 

d 

H2(a) = ^2 + ^ ® h + ia— ® I2 . 

1 uX 
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The operator i is form bounded with respect to —A with relative bound zero. Since V is 
positive, it follows that i I2 is form bounded with respect to H2 with relative bound zero. 

So, BI2(q!) is an entire analytic family in the sense of Kato on D(M2). Furthermore, since H2(q;) is 
unitarily equivalent to IH2 for a G M, we know that ]Hl2(a) is self-adjoint and ct(]HI2) = a(M2{ct)) 
for a G M. Since IHI2 has purely discrete spectrum, we know ]H[2(a) has purely discrete spectrum 
for Q E M. It follows from lemma 4.3 that ]Hl2(a) has purely discrete spectrum V a G C. Since 
112(0) is an entire analytic family in the sense of Kato, the eigenvalues are analytic on C except 
possibly at isolated crossings [19]. Ty(a) unitary implies that the eigenvalues are constant in a 
neighborhood of the real axis and thus crossings will not be an issue. Therefore, the eigenvalues 
are entire functions and constant in a. 

Let P{a) be the projection onto the eigenspace corresponding to the eigenvalue E of El2(a). 
Then P{ct) is entire in a and has the form 



P(a) = — 



\-E\-- 



(H2(a) - A)"^ d\ 



If Q, ao G 



P{a + ao) 



-1 

-1 
27ri 



\~E\ 



(H2(a + ao) - A)"^ d\ 



VF(ao) (E[2(a) - A)"^ W{aQ)-^d\ 



\-E\ = e 

= VF(ao)P(a)VF(ao)"^ . 

For ao G M, the operator valued function /(a) = W{aQ)P{a)W{aQ)~'^ — P{a + ao) is entire 
in a. Since it vanishes V a G M, it is zero V a G C. So P{a + ao) = W{aQ)P{a.)W{aQ)~^ , for 
ao G M and a G C. The hypotheses of O'Connors lemma are satisfied [19]. So, for the eigenvector 
( f\ 

^ = \ , we know ^{a) = W{a)'^ has an analytic continuation to all of C. Therefore 

f, g £ L>(e'^l^l) for any 7 > 0. 

From this it now follows that A/, Ag G L»(eT<^>), for any 7 > 0. To see this, consider 

e7(x> [(^Vn - E) f + Vu g] 
e7(-> [V21 f + {V22 - E) g] 




(^1 



11 



+ 



E) f + e^<"> 1^12 g 

2 

e7(-) V21 f + e^<"> {V22 



E) g 
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Let ;3 > 0. Then, 

J e2-><-) I V21 f\UXdY < II Vi, e-2/3(-> 



,2(x>(7+/3) 



/ 



< 00 



< 00 



So, e^^^Vai/ G L2(m2) and by similar arguments e'^^'-''^ {Vn - E)f, e^^^Vizg, ^'''^ {V22 - E)g 
e L2(r2). Hence, 

■ eT<^') A f \ 
eT<^> A 5 J 

and A/, A5 G L»(e^<^>). 

For V/, Vg G Z)(e'^<^>), we apply Lemma 3.4 of [8]: 

Vp{x) 



Let p G C^{1R'^) and suppose for some C < 00, 

/ ^ + |A/|2)pdx<oo, then 

J ]R 



p{x) 



<2C V X G iR^. If 



\l/2 / /■ \l/2 



+ 



J /I 



We let p (X, y) = e2T<^) . Then 



1/2 



Vp(X, Y) 



\ 1/2 



I A/Mpdx + CM 



1/2 



< 27 V {X, Y) G -K^. We have already shown 



that for / and the right hand side in the lemma is finite for any 7 > 0. So, V/, Vg G D{e^^^'^) 
for any 7 > 0. □ 

Corollary 4.5. Let R{\) = i^2 — for A G p{Ml2)- Let Pe be the projection onto the 
eigenspace associated with E and define r{E) = [{Ml2 — L^)\Ran{i-PE)\~^ > reduced resolvent at 
E. Then, (eT<^> (g) L2) R{X) (g-T^^') (g) L2) and (eT<^> L2)r{E) (e-T<^> ® L2) are hounded on n 
for any 7 > 0. In particular, if ^ G D{e^^'^'^ I2), then R{X) ^, r{E) ^ G D{e^^''^ ® I2). 

Note: See [9] for a proof. 

We need the following lemma for proposition 4.7 

Lemma 4.6. For fixed t G iR", there exist K > A' > and S{t) > 0, such that ifp G IR^ satisfies 
Tr,=iP]>S{t?, then 



< 



kY^p- 



Furthermore, S{t) is uniformly bounded for t in compact subsets of M^' 
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Proof: 

Let 



,2 



1/2 



1/2 



. We prove the Lemma with 



S(t) = l+4||t|| , k = 7/16, and = 17/16. We first show that for this choice of 5 (t), K = 17/16: 

n 



i=i 



1 + 



< 



1 + 



(1 + 4 



Vie / ' 



In particular notice that this argument also shows 



II ii2 II ii2 -^'^ii ii2 

WvW + * < -TT, b 

lb 

Now we show that for this choice of S{t)^ K = 7/16: 



1/2 



> 



> T.(pf-*')-^Y.\p^\\'^ 



3=1 



(4.8) 
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□ 



Proposition 4.7. Let 



e H be a solution of Il2^ = , with E > 0. Then, for 



any 7 > 0, and any a G L>°/, D'^g G D(e^<^>), where = d^' d. 



«2 

X • 



Proof: 

We use a Paley- Wiener Theorem, Theorem IX. 13 of [18]: 
Let (j) G L'^{IR^). Then e^^^ ' </> G L?'{IR^) for all 7 < 7' if and only if has an analytic continuation 
to the set {p : |Imp| < 7'} with the property that for each t G iR" with |t| < 7', + it) G 
L'^{IR^), and for any 7 < 7', sup|^|<^ (?!)(• + it) < 00. 

If a function (/> satisfies the conditions in this theorem we will say that (p is "P-W". Let 
Pj = —i ■ We present the proof for general n. In our case we have n = 2 with xi = X and 
X2 = Y. 

Proposition 4 shows that / and g are P-W for any 7' > 0. In particular we know that 
/, g are analytic everywhere. So the analyticity condition will be a non-issue in the course of 
the proof. V/, Vg, A/, A^f are also P-W for any 7' > 0. So p Pj fip), P ^ PjSip)^ 
p ^ Y^j=i Pj f{p)j ^iid p ^ Yl]=i Pj dip) ^-^^e P-W for all 7' > 0. 

Let S{t) = 1 -|- 4 I |t 1 1 and Bs be a ball of radius S centered at the origin. Since X^j=i pf f{p) 
is P-W, with (4.9) we have 



/ \\pt\f{p + it)?dp < 



sr\B, 



S(t) 



< 00 



itjf 



\f{p + it)\^dp 



(4.10) 



uniformly for t in compact subsets of iR". We only show results involving /. The same results 
hold with / replaced by g. 



Note that since S{t) and /(• + it) 



are uniformly bounded for t in compact subsets of J?", 



we only need to prove estimates for ||p|| > S{t). All of the integral estimates that follow hold 
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uniformly for t in compact subsets of M". From (4.8) and (4.10) we have 

\Pj + \'^\Pk + itk I ^ \f{p + it) P dp 



< 



< 



17 



+ 



16; J]R-\Bs^,^ 



\fip + it)\^dp 



\fip + it)\Up 



< oo . 

It follows that d^^ d^^ f G L»(e^<^>) for any 7 > 0. A gain the same will hold for g. 

We now start an induction on the length \a\ in -D"/ and D"g. Assume that 
D>^f, Dl^g £ D(e'^<^>) for any 7 > and any \(3\ <m-l. It suffices to prove that G L>(e^<^>) 
for any 7 > and any \a \ = m. 

Following the notation in the proof of Proposition 4, the eigenvalue equation gives us 

A/ = Vug + {Vn - E)f , 
^9 = V21 f + {V22 - E)g . 

where Vn, V12 = V21, and V22 are polynomials in xj. Let \a' \ = m — 2. Since the Vij are 
polynomial, our induction hypothesis gives us D"'A/ G i:)(e'^<^>) for any 7 > 0. It follows that 
for jfe G {1,2,--- ,n} 

2 



\Ph + ^tjl I ^ \Pj2 + «*J2 I ^ • • • \Pj^-2 + itjn 



iR"\Bs(t) 

and from (4.9) we have 



|/(p + it)pdp < 00 , 



\Pjl + iijl I ^ |Pi2 + iij2 I ^ • • • \Pjm-2 + i* 



'iR"\Bs(.) 

Since the jk are arbitrary, we have 



Jm-2 



\f{p + it) \^ dp < 00 . 



00 > 



n » 

X] \Ph + I ^ \Ph + it 32 I ^ • • • |Pi™-2 + ^^jm-2 I ^ I |P 1 1^ I /(P + i*) P rfp 
1 ^ in 



Jl J2,--- Jm-2 = 



,1 iiR"\Bs(,) 



juj2,- ,jm-2 = l \^S(t) 



2 ^ 

m-2' 



\f{p + it)\^dp 
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= L i\\p\\' + m\'r-'\\p\\'\fip+^t)\'dp 

JlR\Bsit) 

> [ Hp 11'^"-'^ Hp 1 1' \fiP + it)\^dp 

\fip+it)\^dp 



Then using (4.8), we have for any jk S {1, 2, • • • , n} 



/ \Pjl + itjl I ^ \Pj2 + itj2 I ^ • • • \Pjm + itjm I ^ |/(P + I ^ ^p 

JjR\Bs(t) 



< / „ (IIpII' + \\t\\'r\f{p + ^t)\'dp 

JIR \Bs(t) 

17^ - 



< oo . 

So, for arbitrary G {1, 2, • • • n}, p i— > • • • Pj,n fip) P-W and it follows that D°'f £ 

D{e^^^^) for any 7 > and any \a\ = m. The same argument will work with / replaced by g and 
the proposition is proved. □ 

Lemma 4.8. Let ^ = ^ ^ = (^2 - A)~^ for A G ^(IHIs), and r{E) = (Ma - E)~'^ he 

the reduced resolvent at E. If f, g e and {D" fgi h)^ £ D{e'^^'^'^ (g) I2), for all a G and 
any 7 > 0, then ® /a) i2(A)^', (Z)° ® /2)?^(£^)^' G L'(e'^<^> ® h), for all a G and any 

7 > 0. 

Proof: 

First note that for any 71 > 72 > and j, = 0, 1, 2, • • • , there exists M > such that 

This relative bound implies that if </> G L>(ei'<^>) for all 7 > 0, then Y'' 4> G L'(eT<^>) for all 
7 > 0, and arbitrary j,k = 0, 1, 2, • • • . 

By an argument similar to the one by which we obtained f,g€ C°^{]B?) in the proof of 
Proposition 4.4, R{\) and r{E) map functions from C^iM^) © C^iJR^) to 

The following identity holds as long as the terms on the right hand side are in 
L'^{m?) © L2(iR2): 

{dx ® h) -R(A) $ = R{\) (dx h)^ - R{X) [{dx ® h)(y)] R{X) (4.11) 
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{a + b)X bY 
bY {a-b)X 



where [{dx ® h){V)] = 
compute [dx ® h, RW]'- 

{{dx ® h) R{X) - RiX) (dx ® h)} $ 



{dx ^ h) i ] - R{X) {dx ® h) {m - A) ( 

V ^2 / V ^^2 



To see this, let R{X) ^> 



and we 




R{\) 



R{X) 



dx 








dx 


(;: 










A) dx -01 
A) dx 1^2 



+ 



Viiljji + Vi2 V'2 
V21 V'l + V22 tp2 



dx{Vu)ilJi + dx{Vi2)iJ2 
dx{V2l)Tpl + dx{V22)ll^2 



+ 



Viidx{i'i) + Vi2dx{iJ2) 

V2ldx{'>Pl) + V22dx{llJ2) 




f*' +[(ax««(v)i *' 

dxW2 J \ 1P2 




- R{X) (M2-A) 

R{X) [{dx (E) hW)] 
R{X) [{dx0 I2){V)]R{X)^. 



Clearly (4.11) holds with X replaced by Y. 

From the hypotheses on $ and Corollary 4.5, we know that for all 7 > 0, R{X) {dx (X" I2) ^ £ 
i:>(e^<^> ® I2) C L2(iR2) L'^{]R^). From Corollary 4.5 and the note above, we know that 
R{X) [{dx ^ l2){V)]R{X)'^ G L»(e^<^> I2) C L2(iR2) ^ ^^(iR^) for aU 7 > 0. From this we 
see that (4.11) holds when apphed to ^ and therefore {dx h) R{X) ^ £ D{e''^''^ ® h) for all 
7 > 0. Similarly, {dy ® h) R{X) ^ G L>(e'><^> h) for ah 7 > 0. 

By applying (4.11) repeatedly, we see that (D" (8> I2) R{^) ^ is a linear combination of terms 
of the form 



R{X) [{D'^^ I2){V)] R{X) [{D'^' I2){V)] ■ ■ ■ R{X) [{D' 



l2){V)]R{X) (D°" h)^, 



where YlJLi [oij[ = let |- Since the [(Z?"j /2)(^)] are matrices with polynomial entries, we use 
Corollary 4.5 and the note above to obtain (Z)" (g) I2) R{X) ^ G D{e^^'^'' ® h) for aU q G and 
7 > 0. The conclusion involving (D" I2) r{E) ^ follows by writing the reduced resolvent in the 
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form (see Theorem XII. 5 in [19]) 

r{E) = ^ j R{\) 

ZTTl J\x-E\ = r>0 

□ 



\-E 



Theorem 4.9. For k> 2, let ^('^"2) 



formulas of chapter 3. Then, f^^-'^\ g^^-'^\ Tpv' G C^ilR^), f^^-'^\ ^C'-^) 



(k) 



(k) 

and ip , be determined by the perturbation 



(k) 



an 



d f^^^ 



2) „(fc-2) 



(k) 



el 



(k) 



el 



e L'(e^<^)), for any 7 > 0. In addition, D^/C^-s)^ jja g{k-2) ^ 



G D(eT<^)) for all a G and any 7 > 0. 



Proof: 

We refer to a function in D{e"'^^'^) (or D(e'^^^^ I2) ) for any 7 > 0, as exponentially decaying 
with arbitrary 7. We first note that from the proof of lemma 4.8, multiplication by polynomials 
in X and Y preserves exponential decay with arbitrary 7. 



Since 



.(0) 



is determined at second order in e as an eigenfunction of EI2, we already 



know from propositions 4.4 and 4.7 that ^^^^ satisfies the conclusion. 
The /(^) and 

gW given by equation (3.20) are determined by followed by a projection 
and reduced resolvent EI2 — E^'^\ acting on "i/^^h By corollary 4.5 we know that the reduced 
resolvent preserves exponential decay with arbitrary 7. The projection Q± was the projection 
in 7i onto the subspace perpendicular to the eigenspace of the eigenvalue E^'^^ of EI2. From 
proposition 4.4, we know that the eigenvectors of EI2 have exponential decay with arbitrary 7, 
and so it follows that Q± will preserve exponential decay with arbitrary 7. Since the matrix 
entries of M3 only contain polynomials and derivatives in X and Y, we know from lemma 4.8 that 
(H3 ^f^'^)) will have exponential decay with arbitrary 7. It follows that ^^^^ will have exponential 
decay with arbitrary 7. By a similar argument, ^^^^ G C°° © C°°. From the definition of Q± 
along with proposition 4.7, we see that all of the derivatives of Ha ^'^''^ are exponentially 
decaying with arbitrary 7. It then follows from lemma 4.8 that all of the derivatives of ^'-^^ are 
exponentially decaying with arbitrary 7. 

Recall that = = ip^f^ = 0. From equation (3.21) we know that 



(3) 



< 



+ 



(0) 



(0) 



-1 



P 



(0) 



p 



(0) 



dx 



dx 



(0) 



(0) 





5/(0) 


+ 


el 


dX 




dfio) 


+ 


el 


dX 



ihP±) 



(0) 



p 



(0) f 5£i 
dy 



hP±) 



(0) 



P 



(0) 



dy 



(0) 



(0) 





df(o) 




dY 


el 






dfiO) 


el 


dY 
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By assumption, (^^r ) G Ti-ei, and 
we have 



(0) 



{hP±)^^^ and P^' are bounded operators on Tiei- So 



>(0) 



(3) 



el 



< A 



dX 



+ B 



a/(o) 



dY 



for some positive real numbers A and and 

,,(3) 



(3) 



is exponentially decaying for arbitrary 7 by 



proposition 4.7. Also, G C°°{]R^) since its (X, y) dependence comes strictly from derivatives 



(3) 



el 



is exponentially decaying with 



of /'•'^^ and g^'^K By a similar argument, we see that 
arbitrary 7, from proposition 4.7. 

One can now use induction on k to show the conclusion. For the induction hypothesis, assume 
that and i/'x given by the perturbation formulas of chapter 3 satisfy the conclusions. 

Using equations (3.27) and (3.28) to determine ^(^^"2) and V'j^^) the conclusion follows from the 
propositions and lemmas previously proved. □ 



5 The Eigenstates of the Leading Order Hamiltonian 

We adopt the following notation throughout: 

1. The operator of nuclear angular momentum about the z-axis is denoted by 
.d 



-. The operator of total electronic angular momentum about the z-axis is 



denoted by . The operator of total angular momentum about the z-axis is denoted by 
L™T ^ (L^^c^/) + (I^Lf). 

2. We let L^{x) be the associated Laguerre polynomials, as defined in [15]. 

The first non-vanishing terms in our perturbation expansion are E^'^\ f^'^\X, Y), g^^\X, Y) 
arising from the eigenvalue equation 



i^(2) 





where 



1 



A 



h + 



bXY 



\ 



bXY 



X^ + 



a + b 



Y' 



2 2 

Let {p, <p) be the usual polar coordinates associated with {X, Y). Define the unitary operators 
U, Z : n ^ Hhy: 



U 



cos(0) — sin{(p) 
sin(<^) cos{(j)) 



and 



Z 



1 

7i 



1 1 

i —i 
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Let r = a^l"^ p, 6 = — , and 
a 



Hu = ^ U-^ H2 U 
la 



1 _ ^ 9_ 1 2 
'2 5^~2^9^"^2'' + 



+ 1 



2r2 



/2 + 



r riMC 
l9 -'^z 



Huz 



-1 



19^ 
2 5r2 



(C/Z) 



19 1 . 
2r ar 2 



/2 + 



v 



2r2 
& 2 



2r2 



Both Hu and -fTj/z commute with L\ 



of the form 



,±i|/|. 



,±i|/|. 



'Vi(r) 

' V'2(r-) 



"'^ (8) /2- So, we search for eigenfunctions of these operators 
0, 1, 2, • • • We warn the reader that although / arises here 



as an eigenvalue of L""*^, at this point we should not associate any physical meaning to /. Here 
we are dealing with the operators Hu and Huz, which are related to IH2 by the operations of U 
and Z. The physical meaning of / will become apparent in theorem 6.1. 

We note that ([/~^ EI2 f/)^' = E"^ was the leading order equation obtained by Renner [20], 
which is unitarily equivalent to our leading order equation IHI2 ^ = E^. Renner showed that 
some of the eigenvalues can be solved for exactly, and used regular perturbation theory up to 
second order to approximate the other eigenvalues. These equations have been studied by several 
other authors, for instance [2, 11]. We repeat some of Renner's results here, but we calculate the 
perturbation series to much higher orders, demonstrating that many of the series are diverging 
inside the region of interest. We also illustrate that there is likely a crossing involving the ground 
state eigenvalue of 1H2 near b ~ 0.925a. The ground state appears to be degenerate for < b < 
0.925a and non-degenerate for 0.925a < b < a. 



5.1 The Exactly Solvable I = States 

The / = states (no angular dependence) are exactly solvable. In this case Hu reduces to 



H 



[1=0] 



/ 1 92 I Q 1 + 62 1 

2 9r2 2r dr 2 2r2 



19^ 
2 



1 d 1- 



■ r + 



2r2 / 
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We recognize that the component equations are of the same form as the radial equation for 
angular momentum 1 states of the two dimensional Isotropic Harmonic Oscillator. From the first 
component equation, the eigenvalues and eigenfunctions (non-normalized) are 

EN^ = {2N+ + 2)^Jl + b , I 1, N+ = 0,1,2,- ■■ 

where = (1+6)^/^ r. From the second component equation, the eigenvalues and eigenfunctions 
(non-normalized) are 





= (2iV_ + 2)Vl-6 , I \, =0,1,2,- 

rlLjv_(rl2)e--V2 

where r1 = (1 — b)^^^ r. 

Since IHI2 is unitarily equivalent to y/aHjj, we see these states give rise to eigenvalues and 
eigenfunctions of H2 given by 

En_ = i2N- + 2) Va-b 

(-rlsin((/))L]y_(rl^)e~^-^/^ 
1 , iV-=0,l,2,... (5.1) 
r~- cos((/>) (rl^) e"''- 

where r1 = (a — b)^/^ p, and 

En+ = {2N+ + 2)Va + b 

r+ cos(0)Li (r+2)e-'-+'/2 
\P,^) = { I , iV+ = 0,l,2,... (5.2) 



r+ sin (<?!>) L]^ (r+^) e 



where = {a + 6)^/^ p. 



5.2 The Perturbation Calculation For the 1^0 States 

In this case, Huz reduces to 



o-[±i'i: 



2 dr"^ 2r dr 2 2r'^ 

1 ^2 Id 1 2 (1^1 ±1)^ 

2 dr^ 2r dr 2 2r^ / 



2 I 1 
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Denote the eigenfunctions of H^^^ by I 1 . It is clear that if I | is an 

g[±m^r) ) \ gM^r) 

eigenfunction of H]^^ with eigenvalue then I I ~ I I is an eigenfunction 

of H^j^^^^ with eigenvalue E. So we only need to find the eigenfunctions and eigenvalues of the 

rrim 

We have not been able to solve for the eigenvalues and eigenfunctions in this case exactly. We 
use regular perturbation theory with perturbation parameter b, letting -ffjj'J = -ffg''^ + bV, where 



/ 1 92 Id In (Kl -1) 



rrM 
^0 



2dr^ 2rdr 2 2r2 



18^ _ J_d_ 1 2 + 
2 9r2 2rdr ^ 2'^ ^ 2r^ I 



1 2 / 1 



y = -r , 

One can show using the relative bound found in equation (4.2), that V is relatively bounded with 
respect to i^g'" on H. So, we know that in terms of 6, -ff[['J is an analytic family of type A for 
small b [19]. Therefore, the eigenvalues and eigenfunctions will be analytic functions of & in a 
neighborhood of 6 = 0. 

We expand the eigenvalues and eigenfunctions of in a series in b: 



= E^'\'\ b\ M'^'I'llft) = (5-3) 

fe=0 k=0 



and solve for the coefficients E^'^^\ ^fc^'''' recursively. Here N indexes the energy levels of -ff^/'J 
for fixed |/|. Again from the two-dimensional isotropic oscillator, the eigenfunctions of iifQ^" are 
known exactly. The lowest state is non-degenerate, with eigenvalue and eigenfunction given by 

eT = \1\, <I'I=( ^ |. (5.4) 

The rest of the states are two-fold degenerate, with eigenvalues and eigenfunctions given by 

<S = I <^l" = 2Ar+|/|, 



<:ln = I h iV = l,2,... (5.5) 



rl'l+iLK(r2)e-'^'/2 
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The functions {e^'<^rl'l L^(r2) e'^'/^} 



form a basis for L'^{]R^) by theorem XIII. 64 of [19]. 



K=0,l,2, 



Then for fixed I G Z,, the functions {rl'l L'j^(r^) e~^^ ^^}k=o,i,2,--- form a basis for the projection of 
L?{]B?) onto r-dependent multiples of e*'*^. We can then use the following orthonormal basis for 
the perturbation expansion: 



B 



NM\ 







where the -Bjvjz| are constants of normalization. The matrix elements of the perturbation V in 
this basis can be obtained explicitly [9]. 



5.2.1 The Non-Degenerate Perturbation Calculation 



Recall from (5.4), for fixed |/| 0, the lowest lying eigenvalue of Hq'^' is E^q''^' = |/| (non- 
degenerate). Since H^j^ is an analytic family, we use non-dengenerate, regular perturbation 
theory. Using the Mathematica software package, we easily computed the exact perturbation 
coefficients up to 28th order for the non-degenerate, lowest lying eigenvalue i?*^'''', for several 
values of |/| (See Figure 3). 




Figure 3: A plot of the perturbation series versus b, of the non-degenerate, lowest lying eigenvalue of 
if^'l up to order 28, for |/| = 1, 2, ■ ■ ■ , 8. The dashed curves are (2A^ + 2)\/l - b, for = 0, 1, 2, 3, 
which are the / = states that we have solved for exactly. 



Recall that we are concerned with the case where < & < a, so that < 6 = ^ < 1. 
The functions i5^'^'l'l(6) likely do not exist as eigenvalues of H^^^ if 6 > 1. Seemingly the radii of 
convergence of the series are smaller as |/| increases. It appears that for |/| = 1, 2, and 3 the radius 
of convergence is likely close to 1 (if not larger). For |/| > 4, the series are behaving erratically 
for values of 6 < 1. The |/| = 4 case appears to be borderline, with radius of convergence possibly 
only slightly smaller than 1. This divergent behavior was seen even from the low order coefficients 
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for the larger values of |/|. The singularities are likely caused by avoided crossings between two 
states with the same value of |/|, as suggested in [20, 2, 11]. 

We highlight the crossing between the \l\ = 1 state and the lowest lying / = state near 
b = 0.925. Recall that for I / 0, the eigenvalue of i^l^'J is also an eigenvalue of H^^^^K 

Together these states correspond to a degenerate eigenvalue of the original operator EI2. The 
I = states are all non-degenerate for 6 > 0. So, this crossing implies that the ground state of 
EI2 is degenerate for approximately < 6 < 0.925a and non-degenerate for 0.925a < 6 < 1. 

5.2.2 The Degenerate Perturbation Calculation 

Recall from (5.4), that only the ground state of i/g^" is non-degenerate if |Z| 7^ 0. In the pertur- 
bation calculation described in section 5.2.1, we used regular non-degenerate perturbation theory 
to obtain the perturbation coefficients for these eigenvalues. From (5.5), we have that for fixed 
|/| 7^ 0, I/q'" also has two-fold degenerate eigenvalues of Eq'^^^ = 2N + \l\ for = 1,2, ••• So 
we must use degenerate perturbation theory to calculate the perturbation coefficients of these 
eigenvalues. Recall from (5.5), the degenerate pair of eigenfunctions corresponding to Eq'^''^ are 
given by *2 and ^J^jl. 

Employing degenerate perturbation theory in the usual manner, we find there is splitting that 
occurs at first order (we omit the details). Armed with the proper linear combinations we can 
then proceed as in the non-degenerate case. Using the Mathematica software package, we easily 
computed the exact perturbation coefficients up to 12th order for the first few eigenvalues iiJ^'I'l 
that are degenerate at zeroth order, for several values of |/| (See Figures 4, 5). 



Figure 4: A plot of the perturbation series versus b up to order 12, of the first 8 eigenvalues of H^^^ that 
are degenerate at zeroth order, for |/| = 1. The dashed curves are (2iV + 2)Vl -feand (2iV + 2)Vl + b, 
for = 0, 1, 2, 3, 4, which are the / = states that we have solved for exactly. 



While the splitting is nicely illustrated, we see that all of the series likely have radii of con- 
vergence well below 1. The radius of convergence appears to decrease as |/| or N increase. The 
divergent behavior was seen even at low orders of the perturbation coefficients. 
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Figure 5: A plot of the perturbation series versus b up to order 12, of the first 8 eigenvalues of H^}^ that 
are degenerate at zeroth order, for |/| = 2. The dashed curves are (2iV + 2)Vl-6and (2A^ + 2)\/l + b, 
for = 0, 1, 2, 3, 4, which are the / = states that we have solved for exactly. 

We also used an elementary finite difference scheme to approximate the eigenvalues at several 
values of b, for < b < 1. The results are given in Figure 6. The plot was generated by 
approximating the lowest lying 17 eigenvalues for a fixed b value, then the value of b was changed 
and the lowest 17 eigenvalues were calculated again. This was repeated at steps of A6 = .01 from 
< 6 < .99. Recall that the / = states were exactly solvable. For comparison, the exact values 
of the lowest lying / = states were plotted as dotted curves. We see that the finite difference 
scheme approximates these eigenvalues so well that the dotted curve are hardly distinguishable 
from the finite difference approximation of these eigenvalues. Near 6 = 0, the 17 eigenvalues that 
are being approximated can be identified by their values at 6 = 0: 

1. The curve that has value 1 at 6 = is actually two overlapping eigenvalues of 1HI2 corre- 
sponding to the degenerate pair of lowest lying |/| = 1 states, one for / = 1 and I = —1. 

2. There are three curves that have value 2 at 6 = 0. Two of the curves are the non-degenerate 
/ = states (one increases with b and one decreases with b). The other curve is two 
overlapping eigenvalues corresponding to the degenerate pair of lowest lying \l\ = 2 states, 
one for I = 2 and I = —2. These curves together account for four eigenvalues of IHI2. 

3. There are three curves that have value 3 at 6 = 0. Two of the curves are overlapping 
degenerate \l\ = 1 states, (one degenerate pair increases with b and one degenerate pair 
decreases with b). The other curve an overlapping degenerate pair of lowest lying |Z| = 2 
states. These curves together account for six eigenvalues of IHI2. 

4. There are three curves that have value 4 at 6 = 0. One of the curves is a non-degenerate 
/ = state, one is a degenerate pair of |/| = 2 states, and one is a degenerate pair of |/| = 4 
states. Together these curves account for five eigenvalues of IHI2. 

This plot supports the claim that a crossing occurs involving the ground state near b = 0.925. 
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Figure 6: A plot of lowest 17 eigenvalues of EI2 as a function of on < b < 1, as approximated by 
a finite difference scheme. The dotted curves are 2^1 - b and 2^/l + b, which are the lowest of the 
/ = states that we have solved for exactly. 

While the finite difference scheme is crude, we are incHned to trust the quahtative features of the 
results considering the lowest of the exactly solvable I = eigenvalues were so well approximated, 
even near 6 = 1 as seen in the figure. We note that as b increases from zero, avoided crossings 
involving states with the same value of |^| occur, as well as crossings involving states with different 
values of |/|. When the uppermost curve is involved with such a phenomenon it will appear to 
change behavior suddenly without reason, but this is only because we can only see the lowest 17 
eigenvalues at each b. 

6 Degeneracy of the Quasimode Energies 

The eigenfunctions of HI2 provide the zeroth order states for the quasimode expansion. Recall 



that if 



/(O) 



g 



(0) 



is an eigenfunction of IHI2, we have derived perturbation formulas in chapter 

3 that determine the functions f^'^^X, Y), g^^\X, Y), and ijf\x, Y) that enter in equation 
(3.11) as the asymptotic series 

= 



^i(6X, eY) f^'\x, y)e'= + ^2{eX, eY) J] g^^\x, 6^= + (X, Y)€ 



k=0 



k=0 



k=0 
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= E E (*S'^(^' Y)f^''-^\X, Y) + Y)g^^-^\X, Y)) + (X, Y) , 

fc=0 \ 3=0 j 

oo 

k=0 

where {^i{eX, eY), 'I'2(eX, eY)} is the electronic eigenfunction basis. The f^^^ and g'^'^) have 
no electronic dependence (they are scalar functions) and iIj^^ has both electronic and nuclear 
dependence. 

Recall that for |/| 7^ 0, if I I is an eigenfunction of Hjj)^ with eigenvalue then 



f[-m{r) \ ( (7ll'H(r) 



is an eigenfunction of h\j^^^ with eigenvalue E. So if |/| 7^ 0, we 



have two-fold degenerate eigenfunctions of EI2 of the form (recall r = a}/^p) 




U Z 



1 
71 



^i\l\<t> g[\l\]{j'^ 

^ (gi(|«|-i)<ii /[Kl](r) - e*(l'l+i)''^5[l'll(r)) 



(6.1) 



and 



-iKI</'^[|«|](r) \ ^ / e^^(l'l~i)*/[l'll(r) + e-*(l'l+^)<^5[l'll(r) 

-iKI</' /[l'll(r) / ^ \ -i (e-^(l'l-i)9^ /^l^ll(^^) - e~*(l'l+^)'^5[l'll(r) 



(6.2) 



By taking appropriate linear combinations, these degenerate zeroth order functions lead to two 
orthogonal quasimodes using the perturbation formulas of chapter 3, possibly degenerate (no 
splitting) or non-degenerate (splitting). 

We adopt the following nomenclature: We refer to the eigenfunctions of IHI2 that arise from 
the eigenfunctions of -ff[|'J, where |/| 7^ 0, as -|-|/| states. We refer to the eigenfunctions of EI2 that 
arise from the eigenfunctions oiH\j^^\ where |/| 7^ 0, as states. We refer to the eigenfunctions 
of EI2 that arise from the eigenfunctions of H^^^~^'^ as |/| = states. 

Theorem 6.1. Let L™^ he the operator of total angular momentum around the z-axis and 
< 6 < 1. Then: 
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1. For I 7^ 0, each state generates a quasimode <I>^ of H{e) that satisfies L™"^ = \l\ <^>^. 
The corresponding degenerate state generates a quasimode that satisfies = 
and L^^^ = The and quasimodes are orthogonal, and asymptotic to 
two-fold degenerate eigenf unctions of H{e). We see that linear combinations of the these 
two-fold degenerate ib|/| states also generate valid quasimodes. 

2. Each |Z| = state generates a quasimode that is asymptotic to a non- degenerate eigenfunction 
ofH{e). 

In either case, the zeroth order of the electronic eigenfunction basis vectors ^i(0, 0) and ^'2(0,0) 
are linear combinations of eigenf unctions of Lf with eigenvalues ±1. 

Remark: 

The physical meaning of I is now apparent. It corresponds to the total angular momentum 
about the z-axis of the wave function being approximated. From the proof to follow, it will be 
clear that the zeroth order of a quasimode, can be constructed to satisfy L™-^ $0 = Z™"^ $o- 
In this case it is a linear combination of two states of the form 

r.iTOT_i{fnuc)^+{rel) and %OT+i(f„„c) 

where 



'-'lTOT_l 



1)H 



lTOT_l, 



T nuc ' 
1^, ^ 



[TOTj^l 



TOT 



Proof: 

Since [H[e), L^^'^] = 0, we know that the true eigenfunctions ^'(e) of H(e) can be constructed 
to satisfy L™-^^'(e) = /™^'I'(e), at each e in a neighborhood of 0, for some G Z. This 

implies that L™'^^(^ = -/™^^'(e), since L™^^'(e) = -L™^^(^. We can therefore 
arrange so that the asymptotic series <I>e = Yl'k^o ^^^k satisfies L™-^ <I>e = Z™-^ <I>e at each 
order of e. We then know that each order of the quasimode, and its complex conjugate, are 
eigenfunctions of L™^ with eigenvalues Z™"^ and — respectively. 
We now separate into two cases: 
Case 1: |/| / 

In this case, we have degenerate zeroth order states of the form in (6.1) and (6.2). Regard- 
less of whether splitting occurs, assume that we depart from zeroth order with a correct linear 



combination 



fiO) 



9 



(0) 



a 



+ a 



, so that this leads to a valid quasimode. 



which satisfies L J*^"^ $e = Z™"^ • Then the zeroth order $0 must also be an eigenfunction of 
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LJ*^^ with eigenvalue The <I>o function is given by 

^>o = ^'i(0,0)/(°) + ^2(0,0) 5^°) 



a 



^1(0,0) (^qF(°) + + ^2(0,0) [aG^^^ + /?G(oy 

e^d'l-i)-^/!!'!! (^1(0,0) + i*2(0,0)) + e^d'l+i)'^ ^H'll (^i(0,0) - i*2(0,0)) 

+ /3 e-^d'l-i)-^/!!'!! (^i(0,0) - i*2(0,0)) + e-^(l'l+i)<^5[l'l] (*i(0,0) + i ^2(0,0)) 

We now plug this into the equation L™-^<I>o — /™"^$o = 0, and for \l\ > 2, we project along 

e'(l'l-i)'^ 

tions: 







, and e 


-im+i)<t> g[\i\]^ and obtain the following 


four equa- 


Lf (^1(0,0) + 


^^2(0,0)) = 




-|/| + 1) (^-1(0,0) + 


^^-2(0,0)) 


(6.3) 


Lf (^1(0,0) - 


i ^2(0,0)) = 




-1^1-1) (^1(0,0) - 


^^-2(0,0)) 


(6.4) 


Lf (M/i(0,0) - 


i ^2(0,0)) = 


{l^OT 


+ (*i(0,0) - 


^^-2(0,0)) 


(6.5) 


Lf {^i{0,0) + 


^^-2(0,0)) = 




+ |/| + 1) (^1(0,0) + 


i ^2(0,0)) 


(6.6) 



Equations (6.3) and (6.4) hold as long as a 7^ and equations (6.5) and (6.6) hold as long as 
/3 7^ 0. By combining (6.3) and (6.6) we obtain Z™-^ — |/| + 1 = l'^^^ + \l\ + 1 which contradicts 
our assumption that |/| 7^ 0. So, either a = or /3 = 0. Assume that /3 = and take a = 1, so 
that equations (6.3) and (6.4) still hold. Since L™^$o = ^™^^o, we know that L™^$J^ = 
—l^'^^^Q. Using this equation and projecting along e"*^'''"^)'''' /H'l] and e"**-'''^^^''^ (jfU'l], we obtain 
equations similar to (6.5) and (6.6), but with replaced by — Z™"^: 

Lf (M/i(0, 0) - i ^2(0, 0)) = (-/™^ + |/| - 1) (*i(0, 0) - i ^2(0, 0)) (6.7) 

Lf (^1(0, 0) + i ^2{0, 0)) = (-/™^ + |/| + 1) (^1(0, 0) + i ^2{0, 0)) . (6.8) 

By combining (6.3) and (6.8) we obtain l'^^'^^ = |Z| and these equations now reduce to 

Lf (^1(0,0) - i^2i0,0)) = - (^-1(0,0) - ^^-2(0,0)) (6.9) 

(^1(0,0) + i ^2(0,0)) = ^-1(0,0) + i ^2(0,0). (6.10) 

By repeating the argument with a = 0, /? = 1, we would instead find = — From 

this analysis we see that a = 1, /3 = and a = 0, (3 = 1 are correct linear combinations 
that will generate two orthogonal quasimodes and respectively. These quasimodes satisfy 
l^TOT _ |/|(|)^ and L'^^'^ = —\l\^f and that they are asymptotic to eigenfunctions of 
H{e). We note that $q = $q . Since H{e) commutes with complex conjugation, we have that 
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$^ is also asymptotic to an eigenfunction with the same eigenvalue as <&^. Since quasimodes 
are determined by their zeroth order eigenfunctions through the perturbation formulas of chapter 
3, this imples that $f = $f since = ^q. So, the and $f correspond to a degenerate 
pair and we see that no splitting occurs in the perturbation expansion. As a result, any linear 
combination would be a correct one. The and = $^ generated by the combinations 
a = 1, /3 = and a = 0, (3 = 1 respectively, are the quasimodes that satisfy L™^ = |/| <I>^ 



and L™^$f = -|/|$f. 

If |/| = 1, we take projections of L™^ $o - = along /M, e^^-^^W, and e-^^"^ c/H^H, 



and obtain three equations. By proceeding in a similar manner to the analysis in the |/| > 2 case 
above, we would obtain (^^'^ = lifQ! = l, /3 = and /J*^"^ = — 1 if a = 0, (3=1. In either case, 
we would obtain (6.9) and (6.10) and the desired results follow as in the |/| > 2 case above. 
Case 2: |/| = 

If |/| = 0, we have non-degenerate eigenfunctions of BI2 of the form in equations (5.1) or (5.2). 
In any case, we see that <I>o is real. Then from L^*^^ <5o = ^™^^o and L™^ $0 = -^^^^ ^0, 
it is clear that Z™"^ = in this case. By plugging into L™"^ $0 = and taking projections 
along e**^ and e"*"^ in a manner similar to the |Z| 7^ case, we obtain the same relations for the 
electronic basis vectors at zeroth order given by equations (6.9) and (6.10). □ 

Corollary 6.2. For all b in some interval (0,(5), if e is sufficiently small, the ground state of 
H{e) (corresponding to the R-T pair of states we are considering) is degenerate. 



From our perturbation analysis, we have that for all b in some interval (0,(5), the ground 
state of IHI2 is degenerate, arising from the / = ±1 states of Hjjz- The previous theorem tells us 



then this will correspond to the lowest lying eigenvalue of (e) corresponding to the R-T pair. □ 

Our perturbation calculations suggest that there is a crossing involving this eigenvalue with the 
lowest lying / = eigenvalue, somewhere near b = 0.925. The ground state seemingly corresponds 
to these / = ±1 states for < 6 < 0.925 and corresponds to the non-degenerate, lowest lying 
I = state for 0.925 < 6 < 1. We now prove that the ground state of H2 cannot arise from any 
other |/| states. 

Proposition 6.3. Let < b < 1. Then the ground state of Huz is either an \l\ = 1 state or an 
I = state. 

Proof: 



Proof: 



that these generate a degenerate pair of quasimodes and (^f that satisfy = <I>^, 

(^f = $^ and L™-^$^ = — If the quasimode energy lies below the essential spectrum. 



As previously mentioned, since [Hjj z ., L"^"^^] = we assume ^' 
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Hjjz can be written in the form 

1 / |Z|^T2|/ 



41"' = H'u'z + ^A I • (6.11) 

* Up ±21/ 



We now show that H^^^ must have an eigenvalue below the eigenvalues of -ff^'" if |/| > 2. 
Recall that a (i/^l'l) = a (^^y"), so we only consider i^JJ'J. Since (j/p ±2|/|) > for 

|/| > 2, we know from (6.11) that -ff^l'J > for all |/| > 2. It easily follows that the lowest 

eigenvalue of must lie below the eigenvalues of H^j^^"^ for all |/| > 2. 

The presence of the off-diagonal terms in H^^^ when 7^ 0, implies that both components of 
the eigenvectors must be non-vanishing. Let ^ be the eigenvector corresponding to the lowest 

[21 

lying eigenvalue of Hjj^. Then, 
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We see that H^l}z least one eigenvalue below the eigenvalues of So, the ground state 

of Hjjz must correspond to the ground state of H^^z or h\^^^ . □ 



7 Proof of the Main Theorem 

Here we use the quasimode expansion constructed in section 3 to sketch the proof of theorem 2.1. 
Our candidates for the approximate wave function and energy in the theorem are 

= F{ep) ^,,K, where 

(K-2 K \ 

(*i(eX,ey)/(^)(X,y) + ^>2{eX,eY)g'^^\x,Y)) + ^ e^i;f{X,Y) 

and E^j^ = X^jLo where the f^^\ g^^\ and E^^^ are determined by the pertur- 

bation formulas in section 3. The cut-off function F{ep) is needed to restrict the analysis to a 
neighborhood of the local minimum of the electronic eigenvalues Ei{p) and E2{p) at p = 0, where 
El and E2 are isolated from the rest of the spectrum of h{eX, eY) and also where the functions 
and operators that we have expanded into powers of e (such as ^i{eX, eY)) have asymptotic 
expansions (recall that {x,y) = {eX, eY) and p = ep). We require that the cut-off function 
F{p) : IR^ [0, 1] be smooth in both variables x and y. It has support in some neighborhood 
where p < S. Also, F{p) = 1 for p < i?, where < R < S. So, the derivatives of F{p) with 
respect to x and y vanish outside the region R < p < S. 
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We have 

{H{e) - E,^k) ^e,K = F{ep) {H{e) - E,,k) - y [ ^x,Y , F{ep) ] ^,,k. (7.1) 

To prove the theorem, one can first show the norm of ^e,K is asymptotic to e. It then suffices to 
prove that both terms on the right hand side of (7.1) are finite linear combinations of the form 
e"^ G, where ^ < oo and J > K + 2. Recall from chapter 3 that Pi was the projection 

in Hei onto {^'i, ^'2}^- We can write 

{H{e)-E,^K)^e,K = Xi{e,X,Y)^i{eX,eY) + X2{e,X,Y)^2{eX,eY) + x±{e,X,Y), 

where x± = -Pi [iH{e) — E^ k) ^t,K\ (so xi X2 have no electronic dependence, but xi. does 
have electronic dependence). 

The analysis regarding X2 and x± is similar to that of xi a-iid will be omitted. Using (3.12) 
with our definition of '^e,K-, we have 

2 K-2 K-2 K-2 

j=0 j=0 j=0 

2 4 K-2 4 K-2 

j=0 j=0 

V j=o j=o ^ j j=0 1=0 

To show ||F(e/3) xi i — ^ e^^^i we can consider the terms in the above equation 

separately and use the triangle inequality. Analogous to equations (3.14) and (3.15), we expand 
all functions with (eX, eY) dependence into powers of e, however, we truncate the series here and 
add an error term. For example, we can write hii{eX,eY) = X^jLo ^"'^n "I" hfi 'iX,Y), 
where we know hlY' {X,Y) is in C°°{X,Y) and is bounded by a polynomial in X and Y of order 
+ 1 on supp(-F(ep)). If we do this, we know all terms of order , for j < K, will cancel in the 
above equations, since the terms of /, g, E, and 'tp± were chosen using the perturbation formulas. 
We show how to deal with the hn term arising in F{ep) xi ^1 only, the rest of the terms are 
handled similarly. Considering only expressions of order e-^+^ or higher, this term can be written 
K-2 K-2 

hn{eX,eY) = Yl e^^^ hf} f^^'> + Y ^^^^^^ iX,Y) f^^l 

j=0 0<j<K-2 j=0 

0<1<K 
j+l>K+l 

Then using the results of section 4, in particular that /^-'^ G D{e"'^^'^), one can show that 

K-2 

F{ep)hu{eX,eY) Y f^^'^ ^i{eX,eY) 

j=0 

0<1<K 
j+l>K+l 
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where Cm-, Dm-, Cij, Dij < oo. We see that this term is indeed of order greater or equal to 
Q^^K+2y All of the terms of xi ^^'^ X± can be handled in a similar fashion using the results of 
theorem 4.9. 

The term involving the derivatives of F in equation (7.1) are handled using theorem 4.9 as well. 
The derivatives of F are supported away from the origin and the terms of ^e,K are exponentially 
decaying. We consider the terms involving derivatives with respect to X. Let M be larger than 
supj^2 \ ^^\ and supjj2 Then, using theorem 4.9, one can easily show that 

52 1 

:2 , F{ep) ^e,A' 



< Me 



+ 2e 







dX 


'Hnuc®'Hel / 



< 0(e°°). 

The conclusion of the theorem follows. □ 



Appendix 

We now argue that the odd terms in the E{e) series must be zero. See [9] for a detailed proof 
which utilizes the perturbation formulas derived in chapter 3. 

The Hamiltonian of interest in terms of the scaled nuclear coordinates {X,Y) = {x/e,y/e) is 
given by 

H{e) = ^Ax,y + h{eX,eY), 

where h[x, y) is the electronic hamiltonian that also contains the nuclear repulsion terms. If E{e) 
is an eigenvalue of H{e), then E{—e) is an eigenvalue of 

H{-e) = ^Ax,y + h{-eX,-eY). 
Under the unitary change X = —X, Y = —Y, we see that H{—e) becomes 

H{-e) = ^^xy + h{eX,eY). 

It is clear that H[e) and H{—e) share the same eigenvalues. This does not immediately imply 
E{e) = E{—e), since there could be a pair of eigenvalues related by Ej^{€) = Eb{—()- However, 
this would imply that Ej^^ = E^ , and then theorem 6.1 implies £^A(e) = Esi^)- Therefore, 
E{e) = E{—e), and as a result the odd terms in the expansion must vanish. 



46 



Acknowledgement 

It is a pleasure to thank Professor George A. Hagedorn for his advise and many useful comments. 

This research was supported in part by National Science Foundation Grant DMS-0600944 while 
at Virginia Polytechnic Institute and State University, and also by the Institute for Mathematics 
and its Applications at the University of Minnesota, with funds provided by the National Science 
Foundation. 

References 

[1] M. Born and R. Oppenheimer, Zur Quantentheorie der Molekeln. Ann. Phys. (Leipzig) 84, 
457-484 (1927). 

[2] J. M. Brown and F. J0rgensen, in Advances in Chemical Physics, edited by I. Prigogine and 
S. A. Rice, John Wiley & Sons, Inc., New York, 1983, Vol. 52, p. 117. 

[3] J.-M. Combes and R. Seller, in Quantum Dynamics of Molecules: The New Experimental 
Challenge to Theorists, edited by R. G. Wooley, NATO Advanced Study Institutes Series, 
Series B, Physics v. 57, Plenum Press, New York, 1980, pages 435-482. 

[4] J.-M. Combes, P. Duclos, and R. Seller, in Rigorous Atomic and Molecular Physics, edited 
by G. Velo and A. Wightman, Plenum Press, New York, 1981, pages 185-212. 

[5] G. A. Hagedorn, Ann. Inst. H. Poincare Sect. A. 47, 1-16 (1987). 

[6] G. A. Hagedorn, Commun. Math. Phys., 116, 23-44 (1988). 

[7] G. A. Hagedorn and J. H. Toloza, Int. J. Quantum Chem. 105, 463-477 (2005). 

[8] G. A. Hagedorn and A. Joye, Commun. Math. Phys. 274, 691-715 (2007). 

[9] M. S. Herman, Born- Oppenheimer Corrections Near a Renner-Teller Crossing, Ph. D. Thesis, 
Virginia Polytechnic Institute and State University. 

[10] G. Herzberg, E. Teller, Z. Phys. Chemie B21, 410 (1933). 

[11] G. Herzberg, Electronic Spectra of Polyatomic Molecules (Princeton, New Jersey, D. Van 
Nostrand Company, Inc., 1966). 

[12] P. Jensen, G. Osmann, and P. R. Bunker, in Computational Molecular Spectroscopy, edited 
by P. Jensen and P. R. Bunker, John Wiley & Sons, Inc., New York 2000, p. 485. 

[13] M. Klein, A. Martinez, R. Seller, and X. Wang, Commun. Math. Phys. 143, 607-639 (1992). 

[14] T. J. Lee, D. J. Fox, H. F. Schaefer III, R. M. Pitzer, J. Chem. Phys., 81, 356-361 (1984). 

[15] A. Messiah, Quantum Mechanics (New York, John Wiley & Sons, Inc., 1958). 



47 



[16] M. Peric and S. D. Peyerimhoff, in The Role of Degenerate States in Chemistry: A Special 
Volume of Advances in Chemical Physics, Volume 124, edited by M. Baer and G. D. Billing, 
Series Editors I. Prigogine and S. A. Rice, John Wiley &: Sons, Inc., New York 2002, pages 
583-658. 

[17] M. Reed and B. Simon, Methods of Modern Mathematical Physics I. Functional Analysis 
(New York, London, Academic Press, 1972). 

[18] M. Reed and B. Simon, Methods of Modern Mathematical Physics II. Fourier Analysis, Self- 
Adjointness (New York, London, Academic Press, 1975). 

[19] M. Reed and B. Simon, Methods of Modern Mathematical Physics IV. Analysis of 
Operators (New York, London, Academic Press, 1978). 

[20] R. Renner, Z. Phys. 92, 172 (1934). 

[21] G. Worth and L. Cederbaum, Annu. Rev. Phys. Chem. 55, 127-158 (2004). 
[22] D. Yarkony, Rev. Mod. Phys. 68, 985-1013 (1996). 



48 



